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Abstract 



We study how next-to-leading logarithms modify predictions from 
leading logarithmic soft-gluon resummation in the heavy quark pro- 
duction cross section near threshold. Numerical results are presented 
for top quark production at the Fermilab Tevatron and bottom quark 
production at fixed-target energies. 
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1 Introduction 

The calculation of hadronic cross sections in the elastic limit (i.e., near 
threshold) in perturbative QCD involves contributions from the emission 
of soft gluons. In n-th order QCD one encounters leading logarithmic (LL) 
contributions proportional to (— a"/n!)[ln 2n_1 ((l — z)~ l )/(l — z)] + . There are 
also terms with next-to-leading logarithmic (NLL) contributions proportional 
to (-a n Jn\)^n 2n ~ 2 ({l-zy l )/{l-z)} + . These "plus" distributions are large 
near threshold, z = 1, (the precise definition of z will be given in the next 
section) and resummation techniques were originally developed to sum them 
in Drell-Yan production [I]. As the Drell-Yan process involves electroweak 
interactions it has a rather simple color structure. Gluons are only radiated 
from the incoming quark-antiquark pair in the hadronic collision process so 
the amplitude is a color singlet. 

The first resummation of LL terms (as well as some NLL terms) in the 
heavy (top) quark production cross section was discussed [0, |3| prior to the 
discovery of the top quark at the Fermilab Tevatron 0]. The analysis was 
based on the fact that the LL terms are identical to those in the Drell-Yan 
process and, even though a cutoff was used to define the resummed perturba- 
tion series, this paper pointed out the importance of incorporating resumma- 
tion effects in threshold production of heavy quarks in QCD. Recently other 
groups have applied more sophisticated LL resummation methods to the top 
quark production cross section, c.f. Refs. |5|, |j, which avoid the use of an 
explicit cutoff, and lead to slightly different values for the top quark cross 
section. At present the experimental results cannot discriminate between 
them and there is an ongoing discussion as to which method is theoreti- 
cally superior. The quark-antiquark annihilation channel is the dominant 
partonic process for top quark production in pp collisions at center of mass 
energy y/S =1.8 TeV. However the gluon-gluon channel is more important 
in bottom and charm production near threshold in fixed-target experiments 
with proton and pion beams. The LL resummation method of ]2[] has also 
been applied to 6-quark production at HERA-B |7[] and more generally for 
fixed-target bottom and charm production by hadron beams ||. 

The order a 2 corrections to the Drell-Yan process in || have allowed a 
check of the NLL terms in the resummation formulae, thus the theory is in 
excellent shape. However, this information cannot be used in heavy-quark 
resummation since some of the NLL terms are a consequence of the more 
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complicated color structure in heavy quark production. There is no exact 
calculation of the heavy-quark production cross section at order a A s . Even 
though the NLL terms at order a z s for heavy quark production near thresh- 



old were available in 1990 |]10fl , at the time of the heavy quark resummation 
analysis it was not clear how to resum them. Therefore all work has con- 
centrated on the LL terms. Now this situation has changed. A resummation 
formalism which correctly incorporates NLL resummation near threshold has 



recently been presented by Kidonakis and Sterman [11, 12 1. The authors an 



alytically compared the order-by-order expansion of their results with the 
known NLL terms JTU] and found agreement at threshold in both the quark- 
antiquark channel and the gluon-gluon channel. 

In this paper we will apply this new NLL resummation formalism to 
calculate the top quark production cross section at the Fermilab Tevatron 
and the bottom quark production cross section at fixed target energies. We 
are particularly interested in the size and therefore the phenomenological 
importance of the NLL terms with respect to the LL terms. To determine 
the relative importance of the NLL terms, we will use the cutoff method 
proposed originally in ||. There are several reasons for this, which we now 
discuss. First, as already stated, we are primarily interested in the corrections 
the NLL terms make to the LL terms. Next, by varying the cutoff the failure 
of the perturbative expansion and the onset of the nonperturbative region 
may be studied directly, bypassed in the newer methods. Finally, given 
the complexity of the formulae it is best to use the simplest approach to 
incorporate the NLL terms in a first test of their magnitude. An analysis 
of the NLL terms is very important since the present LL resummations are 
based either on neglecting them entirely |J or retaining only some of them 
|J. These terms are not universal and could possibly be of importance in 
several areas of perturbative QCD besides the analysis of the heavy quark 
cross section, including the production of large transverse energy jets and the 
production of supersymmetric particles. The conclusions of our study should 
not be generalized to these other reactions since each process is different and 
requires a separate study. 
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2 General Formalism 

Here we present some basic formulae which we need for our analysis, based 



on the theory developed by Kidonakis and Sterman [JTT], |i2fl . For heavy 
quark production to be kinematically allowed, the square of the partonic 
center of mass energy, s = x a XbS, must be larger than the threshold value 
of the invariant mass of the heavy quark- antiquark pair, Q 2 = 4m 2 . When 
the heavy quarks are produced with zero velocity, the true threshold and the 
partonic threshold are equivalent. In either case, the plus distributions which 
must be resummed are functions of 

O 2 

' = t- 

Our calculations are based on the factorization of soft gluons from high- 
energy partons in perturbative QCD [T3|, [T4]. All the singular distributions 



in the heavy quark (Q) production cross section can be expressed in the form 
d(J]%\h2 \ ^ \ ^ f dx a dxh 



m2 f —— ( t>a/h 1 { X a,Q 2 )(t>b/h2{Xb,Q' 2 



dQ 2 d cos 6*dy a? ij •* Xa Xb 

xi ( s _l, n |)^»)(_£_,9., a , W2 )), (2. 2) 

where x a , Xb are the momentum fractions of the incoming partons, y is the 
rapidity of the QQ pair, and 9* is the scattering angle in the pair center 
of mass system. The 0's are parton densities, evaluated at the mass fac- 
torization scale Q 2 , in the DIS or MS mass factorization scheme. The hard 
scattering functions a, which depend on the color structure (IJ) of the in- 
teraction as well as on the partonic channel ab and the mass factorization 



scheme, contain all the plus distributions in the threshold region. In |TT 
it was shown that up to NLL it is possible to pick a color basis in which 
moments of the functions a^ b J \z, 6*, a s (Q 2 )) exponentiate with respect to z 
so that 

a { J b J \n,9*,Q 2 ) = C dzz n - l a^\z,9\a s {Q 2 )) 
Jo 

= h T {B* , Q 2 )h*j{6* , Q2) e Eu(n,o*,Q 2 ) . (2.3) 

The hard scattering prefactor is a product of the contributions from the 
amplitude hj and its complex conjugate hj. Therefore / refers to the de- 
composition of the hard scattering amplitude into a color basis, e.g. into an 
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s-channel singlet and octet. The functions hj and hj have no collinear or 
soft divergences at the partonic threshold since these terms have been fac- 
tored into the exponent Ejj in eq. (2.3). The exponential in (2.3), called 



Sjj in [[LI], [12] , satisfies a renormalization group equation with an anomalous 
dimension matrix Tjj. 

The exponents of eq. (2.3) are given by 

1 z n-l _ I 



Eff(n,9\Q 2 ) = - dz— 
Jo I 

x{(2- 



dy Jab)r (/-, v2-r/-i ,\r ,o2\i , (ab) r ^ . w y2 



Z 

h 1 - y 

+gP[a s ((l - ^) 2 g 2 ),r] +^ J) *[a s ((l - z) 2 Q 2 ),9*} } , (2.4) 

where g\, g2 and #3 are functions of the running coupling constant a s with z- 
(and y)-dependent arguments. The parameter r, 1 in the DIS scheme and in 
the MS scheme, effectively changes the lower limit of the y-integral in the MS 
scheme (c.f. jT5j and our later discussion). As it stands, eq. (2.4) is ill-defined 
since it incorporates arbitrarily soft gluon radiation at the point where the 
QCD perturbation expansion diverges. There are different opinions as to how 
to define the exponent in a manner which clearly separates the perturbative 
and non-perturbative regions. We will return to this point shortly. 

To go to NLL in the exponents Eu, we need g\ up to and including 0(a 2 ) 
and g 2 and g% to 0(a s ). Refs. ||, || do not include the NLL g 3 term, treated 
numerically for the first time here. 

The functions g\ and g 2 depend on the mass factorization scheme and the 
identity of the incoming partons but are essentially independent of the color 
structure. The function g\ is scheme dependent because of the r dependence 
of the y integration. It also depends on the identity of the incoming partons 
so that 

(a6)r 1 /~ . „ s 1 at* . 1 _ (ol r 



grWs] = (C a + C b ) I ^ + -K j j , (2.5) 
where C q = Cg = Cf and C g = Ca- The constant K is (c.f. |fL6| ) 



where ri{ is the number of quark flavors. The g 2 contribution is only nonzero 
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in the DIS scheme ||15||, 



3 „ a s 

~0^F — 
I 71 



(2.7) 



Otherwise, 



(qq) MS 



The functions gjp depend on the color structure in the hard scattering but 



gf 



0. 



are mass factorization scheme independent. In [pTT] , |12jl it is shown that they 
are the eigenvalues of the anomalous dimension matrix Tjj which appears in 
the renormalization group equation for the function Sjj. The precise relation 
between the g^ and the eigenvalues of Tjj is 



gjP[a.((l - z) 2 Q 2 ), 6*} = -Aj[a.((l - z) 2 Q 2 ), 6*] 



(2.8) 



showing that gjp and A/ are functions of both the running coupling constant 
a s (evaluated at a ^-dependent scale) and the angle 9* between the directions 
of the incoming parton and outgoing heavy quark. As can be seen below, 



(ab) 



111 



although the eigenvalues are complex, the total g% contribution to E u 
eq. (2.4) is real. 

Consider first the reaction channel q{p a ) + q(Pb) Q(Pi) + Q(Pj) where 
Tjj is two-dimensional. In the s-channel singlet-octet basis the color decom- 
position is into 5 ab 5 ij (singlet) and -5 ab 5 ij /{2N) + 5 aj 5 bi /2 (octet). In this 
basis the components of the matrix Yjj are |12| 



r 2 i 

Tl 2 
^22 



71 

— In (— 

71 11 Ul 
7T C A 11 Ul 




Lp + %7\ 



where 



1 - 2m 2 /s 



In 



VK 



(2.9) 



(2.10) 
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and /3 2 = 1 — Am 2 /s. The Mandelstam invariants for the reaction are s = 
{Pa + Pb) 2 , t = ti + m 2 = (p a - p^ 2 , and u = u x + m 2 = (p b - pi) 2 . 

Tjj is diagonalized in this basis when the parton-parton cm. scattering 
angle 9* = 90° [u\ = ti = —s/2 at threshold) with eigenvalues 

(X 

Ai = A sinE i ct = -C F {L P + 1 + wr) , (2.11) 



TT 



A 9 



A 



— { — Cf(L(3 + 1 + m) + ^ [L, - In (^) + ztt] } 



■m 2 s- 



octet 



It is also diagonalized at partonic threshold s = 4m 2 for arbitrary 9*. 

In the partonic channel g(p a ) + g(Pb) Qipi) + Qi.Pj) the anomalous 
dimension matrix is three dimensional. In the color basis 5 ab 5ji , d abc 7 , " 1 
if abc the components of the matrix Tjj are [TT2 



rii 


Ol s 
TT 


r 2 i 


= o, 


r 3 i 


2a S] 

TT 


r i2 


= o, 


r 2 2 


a s \ 

TT 1 



111 



Ca 
2 



In 



' m 2 s 



ITT 



Tl3 

r 2 3 

T33 



N 2 -4, 
4/V 

2 F31 ' 



T 31 , 



5d 
4 

1^22 



r 



31 



(2.12) 

At 9* = 90° the anomalous dimension matrix becomes diagonal with 
eigenvalues 



Ai 



a 



[C F {L P + 1) + C A m\ 



TT 



^{-C F (L e + l) + §±[L,-ln(^)-»]} 



m 2 s- 



ti 



A 3 = A 2 . (2.13) 
We note that Tu is also diagonal at partonic threshold for arbitrary 9*. 
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3 Heavy-quark production cross sections 

As we have already stated, we are primarily interested in the magnitude of 
the NLL terms, particularly the #3 contribution, relative to the LL results of 
[0, ||, [7|, H . Therefore we use a similar cutoff scheme, modifying the definitions 
of the exponents in eqs. (2.3) and (2.4) to work directly in momentum space 
and avoid working in Mellin space. This method, introduced in 0, exploits 
the correspondence between the moment variable n and logarithmic terms 
in the momentum space variable 1 — z. To directly compare with other 
recent LL resummation schemes in Refs. ||, || would, in principle, require 
three separate studies since each group of authors defines eq. (2.4) differently. 
Therefore, in this paper we only consider the approach based on 0. 
First, we identify 

(3-D 



2m 



where S4 is the invariant mass of the heavy-quark + gluon system which 
recoils against the detected heavy antiquark in inclusive Q production. The 
invariants satisfy s + ti + Ui = S4 so that in the elastic limit S4 — > or 
z — ► 1. Note that eq. (3.1) differs from the identification 1 — z — s^/m 2 in 
Refs. |fj, |5|. The factor of two is necessary to compare the NLL terms JT2 



with the expressions previously given in [JTH] . With the identification in eq. 
(3.1) all the NLL terms in the qq channel are reproduced, not only exactly 
at threshold but also at finite S4 except for one NLL term in C\ which has 
a different coefficient. (This small difference is due to the fact that the QQ 
invariant mass and the Qg invariant mass are different invariants). A similar 
observation holds for the gg channel. 

Further we exploit the fact that a heavy quark- antiquark pair will be pref- 
erentially produced back-to-back in the parton-parton center-of-mass frame 
near threshold. We also choose the case of 6* = 90° to use the results of 
eqs. (2.11) and (2.13). This choice avoids both the diagonalization of the 
three dimensional matrix in the gluon-gluon channel and an extra numerical 
integration over 6*. 

The resummation is done using eqs. (3.17) and (3.19) of [0 where the 
"plus" distributions are removed through integration by parts. Therefore 
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the resummed partonic cross section is defined as 

/ 2^ rs-2ms^ / S4 x da<2(s, Si ,m 2 ) 
a ab {s, m ) = - ds 4 fab r— =■ -, , (3.2) 



Scut 



2m 2 J ds4 



where ah = qq or gg and s cut will be defined below. The other inputs needed 
here are the function f a b, in either the DIS scheme or the MS scheme and the 
differential of the Born cross section. The exponential function f a b{s4 1 /2m 2 ) 
is given by a momentum space version of eq. (2.4). This is best illustrated 
by examining eqs. (3.31) and (3.35) in 0. There we showed how to include 
gluons with an invariant mass s' 4 > s 4 and integrate over their transverse 
momenta up to the kinematic limit. The corresponding result follows by 
first replacing z in eq. (2.4) by z' . Next we replace z' n ~ x — 1 in the numerator 
of eq. (2.4), by —1 and introduce the variables to' — 1 — z' — s' A /(2m 2 ) and 
£ = (1 — y)(l — z')Q 2 /A 2 . By explicitly leaving Q 2 in the following discussion, 
we can later study the scale dependence of the calculation. Note that a s in 
eqs. (2.5)-(2.13) is a function of s'±/2m 2 . 

We first discuss the qq channel. There is only one kinematic structure in 
this channel so that the singlet and octet eigenvalues in eq. (2.11) multiply 
the same function. At this order in perturbation theory however, only the 
octet component contributes, c.f. the discussion in ||11|| . Therefore in the DIS 
scheme at 6* = 90°, the exponential function is 

f%%^)=eME% S + E q - q (\ 2 )), (3.3) 

where 

12 /A2 



2) 



— \ \ — — H oi 2 (e) K 



2 7T ~ S V A2" 



:! '' F «.(^)}- <«) 



Since our calculation is not done in moment space, the u>' integral is cut 
off at ujq = S4/ 2m 2 . Because the running coupling constant diverges when 
u' 2 Q 2 /A 2 ~ 1, the minimum cutoff in eq. (3.2) is s cut = S4 iin i n ~ 2m 2 K/Q. If 
we choose the scale Q 2 = m 2 , s cnt ~ 2mA. In general we choose a larger value 
to stay away from the point of divergence. Once we have chosen a reasonable 
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cutoff, consistent with the sum of the first few terms in the perturbative 
expansion, we can study the size of the NLL terms with respect to the LL 
terms (the order- by-order expansion does not require a cutoff). 
We use the two-loop running coupling constant, 



+ 



Mn(ln£) 



aln£ ' a ln 2 (£) 
UCA-ATfUf 



b = -6 



12tt 

Y1C\ - (6CV + 10CU)?)nf 
(llC A -4T fni y 



(3.5) 



with Ca = 3, Cf = 4/3 and Tj = 1/2. It is obvious that when £ < 1, ln£ 
is negative so that ln(ln£) has a cut and needs a precise definition. We use 
the cutoff on the u/ variable to stop the integration before that point. The 
exact cutoff depends on the quark mass, the scale and A„ f , the QCD scale 
parameter for ri{ light quark flavors, in the parton densities. Both the g\ 
integrals over u' and £ and the g 2 integral over u' as well as the contributions 
to the #3 integration which only contain s' 4 in the running coupling constant 
can be done analytically, as shown in the appendix. 
In the MS scheme, 



(3.6) 



where now 



Note the difference between the upper limits of the ^-integrations introduced 
by the mass factorization scheme in eqs. (3.4) and (3.7). 
The color-dependent g 3 contribution in eq. (2.11) leads to 



e^K) = -[ %W«.fj£-)X = 90°] +KH^f)X = 90°] }. 



<J 2 Q 2 



(3.8) 



in both mass factorization schemes, where i = 1,2. 
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The differential function in eq. (3.2) follows from the kinematic behaviour 
of the Born cross section, represented by 



i$(Mi,«i) = ^ + ^-. (3-9) 



Substituting ti = -l/2{s-s 4 -[(s-s 4 ) 2 -4sm 2 ] 1/2 cos6>*} and u\ = -l/2{s- 
s 4 + [(s — s 4 ) 2 — Asm 2 } 1 / 2 cos 9*}, and taking 9* = 90°, we define, analogous 
to eq. (2.20) in @ 



(o) [(s - s,f - Asm 2 } 1 / 2 B 
F m= 2? m ' 

Differentiating with respect to s 4 , we find 



dF {0) 1 s-8* 

[3(s - s 4 ) 2 - Asm 2 } . (3.11) 



dsi 4s 4 ^( s - S4 )2 _ 4sm 2 
The integrand in eq. (3.2) becomes, for both factorization schemes, 
/ s 4 \ da^js, s 4 ,m 2 ) dF^ 

Hm) d7, = a ql F ~d^~ 

x exp[E qq + E qq {\ 2 )} , (3.12) 

where K qq = 1/N 2 is a color average factor. We note that near threshold 
da^ /ds4 is approximately a factor of two smaller than da^ jds^ in eq. (3.20) 
of Ref. || where the angle 9* was analytically integrated before differentiating 
with respect to s 4 . 

The treatment of the gluon-gluon channel is very similar but now we have 
three distinct color structures. However, only two of them are independent. 
Therefore we define f gg for each eigenvalue so that f gg ^ = exp(E gg + E gg (\j)) 
with a one-to-one correspondence between E gg (\i) and the eigenvalues of eq. 
(2.13) so that 



E gg (x t ) = -£ ^{M« s (^),r = 90°] + A^(^-),r = 90°] } 



(3.13) 
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where i = 1, 2. In the MS scheme, we take 



E, 



mi 



i duo' rQ 2 /^ 2 

w U)' Ju' 2 Q 2 /A 2 £ 



(3.14) 



where K is defined as in eq. (2.6). This choice of K requires some justifica- 
tion, which we have not found explicitly in the literature, due to the fact that 
the discussions of the corrections to the Drell-Yan process are made in the 
DIS scheme. However, it should be possible to prove rigorously that (3.14) 
holds im 



Because there are two kinematical structures, da^/ds^ has two contri- 
butions: 

t± Mi 4m 2 s 

Ml 



ti t\U\ 



1 - 



m 2 s 



F gg,Il( S ^l,Ui) 



l-—^)F ggiI (s,t 1 ,u 1 



(3.15) 



The second kinematical structure, proportional to the first, comes from the 
nonabelian part of the cross section. As in eq. (3.10), we define 

[(a - s 4 ) 2 - Asm 2 } 1 / 2 



-p(o) 



2s 2 "' l ' 

The differential with respect to s 4 of the two components yields 

16m 2 s 



(3.16) 



dFi qli(s,s 4 ,m 2 ) 



1 



+ 



G?S 4 

320mV 



2s 2 y (s — S4) 2 — Asm 2 
1024m 6 s 3 



2(s - s 4 



s 4 



(s - S 4 j J 

dF g g,Il( S i S 4i m2 

dsA 



[s - s 4 

39 



(3.17) 



dF( !!li(s, s 4 ,m 2 ) 



<is 4 



512m 6 s 6(s — s 4 ) 3 64m 4 
2s 2 yf(s - s 4 ) 2 - 4sm 2 |> ~ S4) 3 s 2 s - s 4 

Inserting the remaining color factors yields the gg integrand in eq. (3.2). 



(3.18) 



fgg\ 



2m 2 ' 



ds& 



, S4 = ^ NCFal[[2CFfiii{ j^ (3 , 9) 

^-4)/ 99 , 2 (^)]^ + 4 C ,/ 99 , 2 (^): 



(o) 



^2m 2 



(0) 

Si \ aP gg,II 



dsi 



}. 
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where K gg = (N 2 — 1) 2 is a color average factor. One can numerically 
compare the coefficients of the color terms in da^j /ds± in eq. (3.21) of Ref. 
[0] with the color terms near threshold in da^J/ds4 in eq. (3.19). Those in 
Ref. ]2| are integrated over the angle 9* before differentiation with respect 
to s 4 and are several times larger than those found here. Note that eq. 
(3.19) contains several terms which can have different signs, making it a 
more complicated function of s 4 . Therefore we do not expect our gg results 
to agree with results obtained without attention to the color decomposition 

Because the ui' integral is cut off at s 4 / (2m 2 ), the partonic cross section in 
eq. (3.2) must also have a lower limit of s cut for the integrals to be finite. In 
earlier work up to LL [0, |3|, [7|, || , the lower limit of the partonic cross section, 
so, was taken to be so/m 2 = (/i //i) 2 in the DIS scheme and so/m 2 = (^ /^) 3 
in the MS scheme where A n{ <C fio £t and \i is the renormalization scale. 
It was assumed that the integral over the nonperturbative region below sq 
was negligible compared to the integral over the perturbative region. The 
final resummed cross section was obtained by comparison with the order- 
by-order approximate cross section up to 0(a^) as a function of fi . Since 
each term in the perturbative result is positive, Hq was chosen such that the 
resummed cross section was slightly larger than the sum of the approximate 
cross sections. It was found that for [i = m in the DIS scheme /io ~ (0.05 — 
0.1)m. Larger values were needed in the MS scheme because of the cubic 
power of /zo//^ In the qq channel, /i ~ (0.1 — 0.2)m, while /x ~ 0.35m in the 
gg channel. The needed in the gg channel was larger because of the color 
factor Ca in E gg , see eq. (3.14). As before, we will use the cutoff method since 
increasing s cut away from s 4imin reduces the scale and scheme dependence of 
the resummed perturbative cross section, as we discuss in detail in the next 
section and was also shown previously for the LL resummation [fj, [3| . In the 
next section, we will calculate the fio appropriate for our chosen s cut . 

4 Numerical results 

We first present our results for top quark production. In fig. 1(a) we plot the 
exponents contributing to /^ IS in eq. (3.3) at 9* = 90° as a function of s 4 /2m 2 
with m = 175 GeV/c 2 and ^Js = 351 GeV. Since only the octet component 
contributes to #3 here, only E qg {\2) is shown. We also show, for comparison, 
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the one-loop calculation of g%, E B1S , given in 0. The one-loop and two- loop 
results are quite similar until S4 approaches S4 jm i n . As u'Q/A — > 1, E^ s (gi) 
is larger since E DIS was smoothed near the cutoff ||. The sum E^ s + E qg (X 2 ) 
is always larger than E because the #3 contribution compensates for the 
small, negative g 2 component (|i^ IS (g2)| is shown). The power of 1 — z in a s 
determines the slope of the exponents at small s^/2m 2 . The g 2 contribution, 
linear in 1 — z, is the flattest as s^/2m 2 — > 0. The g\ component shows more 
growth, diverging from E Dls near S4 jm i n , and the #3 contribution, quadratic 
in 1 — z, grows fastest. The growth of the sum is intermediate to that of 
E™ s ( 9l ) and E qq (X 2 ). 

To better illustrate the importance of these contributions to the partonic 
cross section, we show the relative enhancements as a function of s±/2m 2 in 
fig. 1(b). The negative contribution and the slow growth of g 2 near s 4imin 
are reflected in the near threshold behavior of f q q{g\ + g 2 )/ 'fqq(gi) where 
e.g. fqq(gi) = exp(M| IS (gi)) . In the intermediate range of s^/2m 2 , E ms is 
larger than El^ s (gi). Although E qq (X 2 ) dominates /^ IS as s<±/2m 2 — > 1, the 
enhancement in this region is not very significant. Thus the choice of the 
cutoff determines the importance of the NLL contributions. However, a large 
NLL correction is not necessarily reflected in the hadronic cross section since 
it is the parton luminosity relative to the equivalent rj = s/Am 2 — 1 that 
determines the true strength of the corrections. 

Similar results are seen for the exponents in the MS scheme, shown in fig. 
2. The g 2 term vanishes in this scheme. The change in the upper limit of the 
£ integral increases the phase space of E^ s . Both these changes enhance the 
MS result. The g% contribution is identical to that shown in fig. 1. 

In fig. 3 we plot the corresponding results for the gg channel in the MS 
scheme. Note that the real parts of the Aj's are the same in eqs. (2.11) and 
(2.13) so that the strength of the g 3 contributions are the same in the qq 
and gg channels. However, since different color structures in the gg channel 
multiply different derivatives of the Born cross section, there is no single sum 
involving the NLL components as in the qq channel. The E gg (Xi) component 
is quite small and negative near threshold, making the sum E gg + E gg (X 1 ) 
indistinguishable from E gg . The E gg (X 2 ) contribution is the same as in the qq 
channel, resulting in the larger sum E gg + E gg (X 2 ) shown in the dot-dashed 
curve. 

Note that changing the scale effectively shifts the range of s 4 /2m 2 . We 
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found also that, for the same value of S4/2m 2 , the exponents are larger for 
smaller values of Q 2 , as expected. 

We have checked the energy dependence of the g 3 contributions. We find 
that E gg (X 1 ) is always negative, decreasing the NLL correction for all s. Its 
energy dependence is quite strong and thus above threshold it serves to damp 
the resummation. On the other hand, there is very little energy dependence 
of Eu{\2) close to threshold. As s grows large, the \\i{m 2 s / 1\) term in 
eqs. (2.11) and (2.13) dominates the energy dependence, causing Eu(\ 2 ) to 
change sign, eventually leading to a strong damping of the NLL component 
(and probably leading to a smooth energy variation of the resummed cross 
section above threshold). This logarithmic term also has the only explicit s' A 
dependence of the g% contribution at 6* = 90° apart from the argument of 
a s . Since Ai only depends on s' A through a s , these contributions have similar 
behaviours in s. 

In fig. 4 we plot the partonic ti cross section as a function of 77 for the qq 
and gg channels. At i] 1 the gg contribution is numerically smaller than 
the corresponding qq result. The importance of the NLL contributions to 
the partonic cross section depends on the lower limit of the S4 integral. We 
have chosen s cut = 10s4 )in i n , corresponding to fio — 0.26m in the MS scheme 
and 0.13m in the DIS scheme when Q 2 = m 2 , similar to the values used 
previously. Changing s cut changes the effective 77 range. However, note that 
for 77 > 1 , the qq results are nearly scheme independent. The qq result in the 
MS scheme is always larger than the DIS result. We have checked that the 
dependence on s cut is rather weak in this channel. The gg channel is similar 
although the different color structure of the NLL terms leads to a stronger 
variation over all r\ in this channel. 

To illustrate the contribution of the NLL terms to the partonic cross 
section, in fig. 4(b) we show the partonic cross section with only the LL 
contributions for the same value of s cut . At 77 > 1 the qq contributions with 
and without the NLL components are nearly the same but at 77 < 0.02 the 
#3 contribution makes the partonic cross section up to a factor of two larger. 
Both with and without the NLL terms, the scheme dependence increases as 
rj decreases. 

We note that for fixed angle, 9* = 90°, the partonic cross section in the 
gg channel can become negative around 77 ~ 1 — 2. Because dF g °g II /ds4 is 
generally opposite in sign to dF g °g j/ds 4 , when it grows larger than dF g ° g j/ds 4 
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the partonic gg cross section becomes negative. In our calculation of the 
hadronic cross section, we exclude the negative region if it occurs. For the 
particular values of s cut and Q 2 shown here, the partonic cross section is 
always positive. However, for other values of s cut and Q 2 the change in sign 
can occur, both with and without the g^ contribution and thus is not just a 
feature of the NLL contribution at this angle. 

The results in fig. 4 imply that the relative importance of the NLL con- 
tributions to the hadronic cross section depend on which regions in rj are 
weighted most heavily by the convolution with the parton densities. With 
equivalent parton luminosities, if rj ~ 0.5, then we expect the gg contribution 
to be comparable to that of the qq channel while if r\ « 2 then we expect the 
qq channel to dominate. Thus the parton luminosity in the x space probed by 
the hadronic cross section determines the ultimate importance of the NLL 
contribution. A comparison with previous results obtained by integrating 
over the angle 9* |2|, |3| is inconclusive since here we work at fixed 9*. Finally 
we note that previously [§, [| 0, H different cutoffs were chosen in the qq 
and gg channels to match the order-by-order contributions, as described at 
the end of the previous section, whereas here we keep the same s cut for both 
channels. 

The hadronic cross section calculated to NLL is 



^NLL, 



A A J TO J T X X 

where aij{rS,m 2 ) is the partonic cross section, eq. (3.2) and r = (m + 
\/m 2 + s cut ) 2 / S. We evaluate the parton densities at jj? = Q 2 . Since the 
parton densities are only available at fixed order, the application to a re- 
summed cross section introduces some uncertainty. We have used the MRS 
D-' DIS densities for the qq DIS channel and the MRS D-' MS densities 
[ |I~8| , |19|| for both channels in this scheme^. Note that A5 = 0.1559 GeV for 



both sets. The result is not strongly dependent on the parton distributions 
in the x range probed at this energy, x ~ 0.2 at y = for m = 175 GeV/c 2 , 
entering primarily through the value of A5 since different sets of parton den- 
sities have a different A 5 , changing S4 iin i n as well as the value of a s . Taking a s 
to only one loop, as in H, increases ot s and thus the hadronic cross section. 



*The x dependence of the parton densities in the two schemes is very similar. 
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For a more complete discussion of the parton density and a s dependence, see 

In fig. 5(a) we plot the it production cross section at 8* = 90° as a function 
of top quark mass at the Fermilab Tevatron with v^S* = 1.8 TeV in the DIS 
and MS schemes for the qq channel and the MS scheme for the gg channel 
with s cut = 10s 4imin and Q 2 = m 2 . The total MS cross section at 9* = 90° 
is 5.3 pb at m = 175 GeV/c 2 , similar to the angle-integrated results of Ref. 
0. Previous NLO results showed that for m = 100 GeV/c 2 , even though 
?7max — 80 the total hadronic cross section was obtained already for rj ~ 3 
(c.f. fig.l in [RJ). When m = 175 GeV/c 2 , ?7 max = 26. Comparing figs. 4a and 
5a, we can infer that in this case the important rj region is rj »s 0.3. At this 
7] the qq luminosity is six times larger than the gg luminosity. 

To quantify the enhancement in the it cross section at fixed angle due to 
the NLL terms, in fig. 5(b) we show our results with the same values of s cut 
and Q 2 but without the NLL contributions, i.e. no g% component. The NLL 
contribution increases with top mass at this energy. The DIS cross section 
is enhanced between 25 and 30% as the top mass increases from 140 to 200 
GeV/c 2 . The enhancement is somewhat larger for the MS scheme, between 
36 and 43% in the qq channel and 47 to 55% in the gg channel. The change 
in the amount of enhancement is, in part, due to the energy dependence of 
the #3 contribution. For our choice of s cut , the it cross section at m = 175 
GeV/c 2 without the NLL g 3 contribution is 3.7 pb. Thus the NLL terms 
enhance the total MS cross section at 6* = 90° by 43% at this mass and 
value of s cut in our calculation. The relatively large contribution from the 
NLL terms verifies that the lower t] region is most important. 

We have further investigated the s cut and scale dependence of our cal- 
culation. In Tables 1-3 we give the numerical values of the hadronic top 
production cross sections at LL and NLL with s cut = 10s4 )m i n , 5s4 )m i n and 
20s 4imin respectively. For each value of s cut we show the results for Q 2 = m 2 , 
4m 2 and m 2 /4. In the tables we have also indicated the equivalent value of 
the cutoff /i defined in [@ in both the DIS and MS schemes. Note that in 
Tables 1 and 3, the scale dependence of the gg channel is weak compared to 
the qq contributions. When s cut = 10s4 )min , the gg cross section changes 41% 
between Q 2 = m 2 /4 and 4m 2 at m = 175 GeV/c 2 . In contrast, the qq MS 
cross section changes 47% and the DIS cross section by 65%. At the same 
mass but with s cut = 20s4 im i n , the gg cross section changes 10% between the 
highest and lowest scales while the qq MS cross section changes 35% and the 
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DIS cross section 52%. On the other hand, for the lowest s cu t the gg cross 
section has the largest scale dependence, more than a factor of two change 
while the MS cross section changes 53% in the qq channel. 

Note that for a fixed value of s cut , the gg cross section is actually largest 
when Q 2 = 4m 2 than at lower Q 2 . This result is counterintuitive since one 
generally expects that increasing the scale decreases the cross section. This 
is true in the qq channel and also in the gg channel when rj > 1. At fixed 
s 4 , the exponents are larger for lower values of the scale due to the running 
of a s . However, the difference in the factors f gg> i and f gg ^ at low r\ shifts 

the relative weights of dF j/ds^ and dF gg II /ds$ so that the partonic cross 
section is largest for the largest scale studied, Q 2 = 4m 2 , when 77 < 0.5. In 
addition, changing the scale for fixed s cut /s4 im m effectively changes s cu t in 
eq. (3.2) since s 4jmin = 2m 2 K/Q. Therefore increasing Q 2 correspondingly 
increases the effective s 4 integration region. This latter effect seems to be 
the most important because the largest scale also produces the largest gg 
contribution to the cross section when only the LL terms are considered. 

When s cut = 5s4 jiri i n , the gluon contribution is enhanced and actually 
dominates for Q 2 = m 2 and 4m 2 , seen in Table 2. This is not surprising 
since at s cu t = 5s 4imin , /i ~ 0.1m in the DIS scheme and 0.2m in the MS 
scheme when Q 2 = m 2 . At such low values of /j,q in the gg channel, the gg 
cross section can become large, see e.g. figs. 12-14 in 0. Additionally, as 
s cut is reduced, all the contributions to the hadronic cross section increase 
because as s cu t approaches S4 min the exponents grow more rapidly. The 
fastest growth of the exponents in this region occurs when Q 2 = m 2 /4. 

A comparison of the NLL and LL results in the tables helps to clarify 
where the NLL enhancements are largest. As discussed above, the enhance- 
ment for a fixed s cu t is largest for Q 2 = 4m 2 since the larger scale increases 
the effective S4 probed at larger Q 2 . The enhancement is also larger for 
lower Scut due to the running of the coupling constant. Finally we note that 
the NLL enhancement increases with mass since it production is closer to 
threshold at fixed energy with larger masses. 

We conclude our discussion of top production by noting that the scheme 
dependence is rather strong, especially for lower s cu t- The scheme depen- 
dence arises because f^ s > /^ IS , particularly figs. 1 and 2 
and also e.g. figs. 12 and 13 in Ref. [Q]. For s cut = 10s 4jmin and Q 2 = m 2 , 
the MS cross section is 62-71% larger than the DIS cross section at NLL and 
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49-55% larger at LL. The scheme dependence is largest for small s cut and 
large Q 2 because smaller values of s cut are increasingly sensitive to the upper 
limit of the £ integral in E ab (gi), see the appendix. The scheme dependence 
increases slightly with mass at fixed energy. Thus close to threshold the 
scheme dependence is unavoidable. Even though the partonic cross section 
is only weakly scheme dependent for i] > 1, as seen in fig. 4, some scheme 
dependence will remain until new sets of parton densities which have incor- 
porated resummation effects before being fitted to data are available. In the 
absence of such densities we favor the MS results as they have a more reliable 
theoretical basis and the qq and gg channels can be treated consistently only 
in this scheme. 

Since a s depends strongly on the heavy quark mass scale, we repeat 
our calculations for bottom quark production with m = 4.75 GeV/c 2 at 
\/s = 9.51 GeV. The smaller quark mass and the larger A4, 0.23 GeV for 
the MRS D— ' distributions, result in a reduced range in s^jlvr? before the 
exponents diverge, as shown in figs. 6-8. The exponents tend to be somewhat 
smaller for the bottom quark except near S4 >m i n due to the faster running of 
the coupling constant. The running of a s also produces faster growth of the 
exponents at lower Q 2 . Thus the scale dependence is also stronger for bb 
production. We note that the energy dependence of g 3 is stronger for the 
lighter quark mass. While E gg (Xi) is smaller near threshold, it increases 
faster with energy than at the top quark mass. In general, the bottom quark 
cross section is more sensitive to s cut and Q 2 than the much more massive 
top quarks, see also 0, |J. 

The partonic cross sections with s cut = 1.4s 4)min and Q 2 = m 2 are shown 
with and without the NLL #3 terms in fig. 9. This value of s cu t corresponds 
to /io ~ 0.37m and 0.51m in the DIS and MS schemes, somewhat larger than 
those in 0, f|- Varying s cut we find that for 77 < 1, the partonic cross section 
can change by an order of magnitude near 7] ~ 0.1 and a factor of two to three 
at i] ~ 0.5. The largest variation is in the gg channel. Of course the parton 
luminosity in the x range probed at a given energy determines the sensitivity 
of the hadronic cross section. At higher energies, away from threshold, the 
sensitivity to s cut decreases. Similar trends can be seen in charm production 
although for cc production, S4 >m i n already corresponds to the /io values used 
in ||, suggesting that a stable cross section is even more difficult to obtain 
for charm production. 

In fig. 10 we show the bb production cross section at 9* = 90° as a function 
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of beam momentum for pp interactions with s cut ~ lAsj ^j , and Q 2 = m 2 , 
with and without the NLL contributions. As expected, gg fusion is dominant 
|0, [ij. These energies, 20 < y/S < 45 GeV, correspond to 3.4 < r/ max < 21.4. 
At the HERA-B energy, y/S = 39.2 GeV, we find a total MS cross section 
at 6* = 90° of 8.6 nb at LL and 19.1 nb at NLL, more than a factor of two 
enhancement. Note that the bb cross section is somewhat smaller than found 
previously [0, §] , in part because we are only considering a fixed angle. 

The sensitivity to s cut and the scale is detailed in Tables 4-6 for s cut = 1.2, 
1.4, and l-6s4 )TT1 i n and Q 2 = m 2 , 4m 2 and m 2 /4. The same trends can be ob- 
served as we remarked upon for ti production although the enhancements 
are generally larger for bb production. We note that for the values of s cnt 
chosen here, the NLL enhancement and scheme dependence decrease with 
increasing energy since the higher energy probes the region where the correc- 
tions are less important. There seems to be a general decrease of the scheme 
dependence with energy for the bottom quarks. However, it is difficult to 
fairly compare the scheme dependence of top and bottom production since 
our values of s cut do not correspond to the same fio/m ratio in both cases. 
Our previous work on LL resummation favored a fixed ratio to describe pro- 
duction of all heavy quarks J?], || . Thus, while our bottom quark results are 
illustrative only, they clearly indicate that a full NLL calculation including 
all angles is needed to clarify bottom quark resummation. 

5 Conclusions 

We have investigated the numerical importance of the NLL terms in the re- 
summation of subleading soft gluon contributions near threshold for heavy 
quark production. We have shown that, for the resummation method of [0 
and at 6* = 90° these contributions are either numerically small or there 
are cancellations between them. Therefore the inclusion of the NLL terms 
leads to only a moderate change in the top quark production cross section 
at 9* = 90°. Bottom production is much more sensitive to s cu t, making 
definitive statements about the resummed cross sections more difficult, par- 
ticularly in a model without a calculable perturbative cutoff. 

We would like to thank G. Sterman and E. Laenen for helpful discussions. 
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A Appendix 

Here we present the analytical results for the g±, #2 and g 3 (where appropri- 
ate) integrals. As previously discussed, we use the two-loop running coupling 
constant, 



= 



1 
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with Ca — 3, C F — 4/3 and Tf = 1/2. It is obvious that when £ < 1, ln£ is 
negative so that ln(ln£) has a cut and needs a precise definition. We use the 
cutoff on the uj' variable to stop the integration at that point. 
In general, 
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2Ca/7t and = 3Cf/2tt for gg in the DIS scheme and otherwise. We will 
give the results for single integrals over u/ first. For convenience, we define 
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Note that the contributions dependent on L2 set the lower limit on s cut since 
the value of u where u ~ A/Q is the lowest uo at which the integrals 
diverge. 
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Then, in the DIS scheme, the g 2 result is 



3Cf 
2 na 



lnL - — (lnL + 1) - InLi + — (I11L1 + 1) . (A.4) 
L Li 



For the terms in \ which only contain an s' A dependence in the running 
coupling constant, the g 3 integral can also be done analytically. In these 
cases, 



Eab(gs) oc 7^[lnL - -^(lnL + 1) - lnL 2 + -^(lnL 2 + 1) 
la 1 L L 2 



(A.5) 



Thus the g% contributions diverge at a larger value of uq than do the g 2 . 
After the £ integral, the g\ exponent is 
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In the DIS scheme the final result is 
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The result is somewhat simpler in the MS scheme since does not depend 
on uj' . In this case, 
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-^L (ln L - 1) + - In L (ln L + 2) 

+ iL 2 (lnL 2 -l)-^lnL 2 (lnL 2 + 2) (A.8) 

In L + -^-(2 In L + 3) + (18 In 2 L + 30 lnL + 19) 

+\ lnL 2 - -^-(2 In L 2 + 3) - (18 In 2 L 2 + 30 lnL 2 + 19)}" . 



25 



References 

G. Sterman, Nucl. Phys. B281, 310 (1987); S. Catani and L. Trentadue, 
Nucl. Phys. B327, 323 (1989). 

E. Laenen, J. Smith, and W.L. van Neerven, Nucl. Phys. B369, 543 
(1992). 

E. Laenen, J. Smith, and W.L. van Neerven, Phys. Lett. B321, 254 
(1994). 

F. Abe et ai, CDF Collaboration., Phys. Rev. Lett. 74, 2626 (1995); S. 
Abachi et al, DO Collaboration, Phys. Rev. Lett. 74, 2632 (1995). 

E.L. Berger and H. Contopanagos, Phys. Lett. B361, 115 (1995); Phys. 
Rev. D54, 3085 (1996); ANL-HEP-CP-96-51, |hep-ph/960642Tl . 

S. Catani, M.L. Mangano, P. Nason, and L. Trentadue, Nucl. Phys. 
B478, 273 (1996); Phys. Lett. B378, 329 (1996). 

N. Kidonakis and J. Smith, Mod. Phys. Lett. All, 587 (1996). 

J. Smith and R. Vogt, LBL-38282, ITP-SB-95-60, |hep-ph/96093"HB] to 
be published in Z. Phys. C. 

W.L. van Neerven and E.B. Zijlstra, Nucl. Phys. B382, 11 (1992); T. 
Matsuura, R. Hamburg and W.L. van Neerven, Nucl. Phys. B359, 343 
(1991). 

R. Meng, G.A. Schuler, J. Smith and W.L. van Neerven, Nucl. Phys. 
B339, 325 (1990). 

N. Kidonakis and G. Sterman, Phys. Lett. B387, 867 (1996); contribu- 
tion presented at the Les Rencontres de Physique de la Vallee d'Aoste, 
La Thuile, Aosta Valley, March 4-9, 1996, |iep-ph/9607222l . 



N. Kidonakis, SUNY at Stony Brook Ph.D. Thesis (1996), |ep= 



ph/9606474l , UMI-96-36628-mc (microfiche), 1996; N. Kidonakis and G. 
Sterman, EDINBURGH 97/3, ITP-SB-97-24, |hep-ph/ 9705234] , 1997. 



26 



[13] J.C. Collins, D.E. Soper and G. Sterman, in Perturbative Quantum 
Chrodynamics, ed. A.H. Mueller (World Scientific, Singapore, 1989), p. 
1. 

[14] J.G.M. Gatheral, Phys. Lett. B133, 90 (1983); J.C. Collins, in Per- 
turbative Quantum Chrodynamics, ed. A.H. Mueller (World Scientific, 
Singapore, 1989), p. 573. 

[15] H. Contopanagos, E. Laenen and G. Sterman, Nucl. Phys. B484, 303 
(1997). 

[16] J. Kodaira and L. Trentadue, Phys. Lett. B112, 66 (1982). 
[17] G. Sterman, private communication. 

[18] A.D. Martin, R.G. Roberts, and W.J. Stirling, Phys. Lett. B306, 145 
(1993). 

[19] H. Plothow-Besch, 'PDFLIB: Nucleon, Pion and Photon Parton Density 
Functions and a s Calculations', Users's Manual - Version 4.16, W5051 
PDFLIB, 1994.01.11, CERN-PPE; Comp. Phys. Comm. 75, 396 (1993). 



27 



Figure Captions 

Figure 1. (a) We show the contributions to /^ IS in top quark production 
for m = 175 GeV/c 2 and ^Js = 351 GeV as a function of s 4 /2m 2 . The 
solid curve shows ln(/^ IS ) for the octet component at 90°, the dot-dashed 
curve, E q g(X 2 ), the dashed curve, E^ s (gi), and the dot-dot-dashed curve, 
\E q ^ s ( y g 2 )\. The dotted curve shows E DIS as defined in eq. (3.30) of Ref. 
[[J. (b) The ratios of the exponential factors are given to show the enhance- 
ment relative to f qq -(gi) where fqq(gi) = exp(£'^ IS (5f 1 )). The solid curve 
shows fg^ S /fqq(gi), the dashed curve, f qq -(gi + g2)/fqq(gi), and the dotted 
curve, fggigx + g 3 (X 2 )) / fqq(gi) ■ The ratio f qq -(gi)/ exp(£ DIS ), given in the 
dot-dashed curve, shows the enhancement of g\ from eq. (2.5) relative to the 
one-loop value in [||. 

Figure 2. (a) We show the contributions to in top quark production for 
m = 175 GeV/c 2 and -^i = 351 GeV as a function of s^/2171 2 . The solid curve 
shows ln(/^ s ) for the octet component at 90°, the dot-dashed curve, E q q(X 2 ), 

and the dashed curve, E^ s (gi). The dotted curve shows E MS as defined in 
eq. (3.34) of Ref. 0. (b) The ratios of the exponential factors are given to 
show the enhancements. The solid curve shows f^ s /fqq(gi) while the dashed 
curve is the ratio f q q(gi)/ exp(E MS ) where f q q{gx) = expiE^igJ). 

Figure 3. (a) We show the contributions to f gg in top quark production for 
m = 175 GeV/c 2 and y/s = 351 GeV as a function of s 4 /2m 2 . The solid 
curve shows E gg ; the dashed and dotted curves give the contributions from 
the eigenvalues, E gg {\2) and \E gg (Xi)\ respectively. The dot-dashed curve 
shows the sum E gg + E gg (\2); the sum E gg + E gg (Xi) is indistinguishable 
from the solid curve. The dot-dot-dashed curve shows E MS for the gluon 
from Ref. 0. (b) The ratios of the exponential factors are given to show 
the enhancements relative to f gg (gi) = exp(E gg ). The solid curve shows 
fgg{g\ + gsi^i))/ fgg(gi) and the dashed curve, f gg (g 1 + g?,^)) / f gg (gi)- The 
ratio of f gg (gi) to exp(i? MS ), given in the dot-dashed curve, compares the g\ 
from eq. (2.5) relative to the one-loop value in j2|. 
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Figure 4. The partonic top quark production cross section is shown as a func- 
tion of 1] = (s — Am 2 ) /(Am 2 ) for the qq channel in the DIS (solid) and MS 
(dashed) schemes and the gg (dot-dashed) channel. We use s cut ~ 10s4 )in i n 
with (a) and without (b) the NLL g 3 contributions respectively. 

Figure 5. The hadronic it production cross section at 9* = 90° is given as a 
function of top quark mass for pp collisions at the Tevatron energy, v^S* = 1.8 
TeV. We use the MRS D— ' parton densities in the DIS scheme for qq an- 
nihilation (dashed) and in the MS scheme for the qq (dot-dashed) and gg 
(dotted) channels. The sum of the qq and gg channels in the MS scheme is 
given in the solid curve. The results are given for s cut w 10s 4imin with the 
NLL g 3 contributions in (a) and with these excluded in (b). 

Figure 6. We show the contributions to /^ IS and the enhancement factors 
in bottom quark production for m = 4.75 GeV/c 2 and ^/s = 9.51 GeV. The 
labels are as in fig. 1. 

Figure 7. We show the contributions to and the enhancement factors 
in bottom quark production for m = 4.75 GeV/c 2 and y/s = 9.51 GeV. The 
labels are as in fig. 2. 

Figure 8. We show the contributions to f gg and the enhancement factors in 
bottom quark production for m = 4.75 GeV/c 2 and y/s = 9.51 GeV. The 
labels are as in fig. 3. 

Figure 9. The partonic bottom quark production cross section is shown as 
a function of rj. The labels are as in fig. 4. We use s cut ~ 1.4s 4imin with the 
NLL g 3 contributions in (a) and without them in (b). 

Figure 10. The hadronic bb production cross section at 9* = 90° is given as 
a function of beam momentum for pp interactions. We use the MRS D— ' 
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parton densities in the DIS scheme for qq annihilation (dashed) and in the 
MS scheme for the qq (dot-dashed) and gg (dotted) channels. The MS sum 
is given in the solid curve. The results are given for s cut ~ 1.4s 4imin with the 
NLL g 3 corrections in (a) and without them in (b). 
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NLL 


LL 


m (GeV/c 2 ) 


^ IS (Pb) 


^ S (Pb) 


°99 (Pb) 


^ IS (Pb) 


< S (Pb) 


°99 (Pb) 




Q 2 = m 2 = 0.13m, /i^ s = 0.26m 


150 


4.91 


8.04 


3.99 


3.91 


5.85 


2.68 


155 


4.21 


6.94 


3.21 


3.32 


5.01 


2.15 


160 


3.63 


5.99 


2.61 


2.87 


4.33 


1.75 


165 


3.13 


5.20 


2.12 


2.47 


3.74 


1.40 


170 


2.71 


4.51 


1.73 


2.12 


3.23 


1.14 


175 


2.34 


3.92 


1.41 


1.84 


2.80 


0.93 


180 


2.03 


3.41 


1.15 


1.59 


2.43 


0.75 


185 


1.77 


2.98 


0.94 


1.37 


2.11 


0.62 


190 


1.53 


2.59 


0.77 


1.19 


1.84 


0.50 


195 


1.33 


2.27 


0.63 


1.03 


1.60 


0.41 


200 


1.16 


1.99 


0.52 


0.90 


1.39 


0.34 




Q 2 = 4m 2 _^ ^ms = o.igm, = .4i m 


150 


3.83 


6.60 


4.94 


3.06 


4.74 


3.11 


155 


3.30 


5.71 


3.95 


2.61 


4.07 


2.49 


160 


2.83 


4.92 


3.22 


2.24 


3.50 


2.01 


165 


2.44 


4.27 


2.61 


1.93 


3.03 


1.62 


170 


2.10 


3.69 


2.13 


1.65 


2.61 


1.32 


175 


1.82 


3.21 


1.74 


1.43 


2.26 


1.08 


180 


1.56 


2.78 


1.42 


1.23 


1.95 


0.88 


185 


1.36 


2.42 


1.16 


1.06 


1.69 


0.71 


190 


1.18 


2.11 


0.96 


0.92 


1.47 


0.59 


195 


1.03 


1.84 


0.79 


0.80 


1.28 


0.48 


200 


0.89 


1.61 


0.65 


0.69 


1.11 


0.39 




Q 2 = m 2 /4 fi™ = 0.094m, /i^ s = 0.16m 


150 


6.20 


9.53 


3.51 


4.97 


7.13 


2.51 


155 


5.36 


8.28 


2.84 


4.26 


6.12 


2.03 


160 


4.62 


7.17 


2.29 


3.67 


5.29 


1.62 


165 


4.00 


6.24 


1.85 


3.16 


4.57 


1.31 


170 


3.46 


5.41 


1.51 


2.73 


3.97 


1.06 


175 


3.00 


4.72 


1.23 


2.36 


3.44 


0.86 


180 


2.61 


4.11 


1.00 


2.04 


2.98 


0.70 


185 


2.27 


3.59 


0.82 


1.77 


2.60 


0.57 


190 


1.99 


3.16 


0.68 


1.54 


2.26 


0.47 


195 


1.73 


2.75 


0.55 


1.34 


1.97 


0.38 


200 


1.51 


2.42 


0.46 


1.16 


1.72 


0.31 



Table 1: The hadronic tt production cross sections in pp collisions at v^S* = 
1.8 TeV for s cut = 10s 4imin and 9* = 90°. The effective /x from M is given 
for each Q 2 . 
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NLL 


LL 


m (GeV/c 2 ) 


^ IS (Pb) 


°JP (Pb) 


°99 (Pb) 


^ IS (Pb) 


< S (Pb) 


°99 (Pb) 




Q 2 = m 2 ^is = o.094m, /ijf 3 = 0.21m 


150 


5.43 


10.81 


17.31 


4.13 


7.07 


9.67 


155 


4.65 


9.29 


14.13 


3.55 


6.09 


7.89 


160 


4.01 


8.08 


11.66 


3.05 


5.25 


6.44 


165 


3.47 


7.02 


9.54 


2.63 


4.55 


5.27 


170 


2.99 


6.10 


7.93 


2.25 


3.93 


4.30 


175 


2.59 


5.32 


6.53 


1.95 


3.41 


3.54 


180 


2.25 


4.63 


5.38 


1.69 


2.96 


2.90 


185 


1.95 


4.04 


4.47 


1.46 


2.58 


2.41 


190 


1.70 


3.56 


3.68 


1.26 


2.24 


1.99 


195 


1.48 


3.09 


3.06 


1.10 


1.96 


1.65 


200 


1.28 


2.71 


2.55 


0.96 


1.71 


1.37 




Q 2 = 4m 2 /iP IS = 0.13m, /i^ s = 0.33m 


150 


4.15 


8.80 


26.50 


3.20 


5.68 


14.01 


155 


3.55 


7.58 


21.69 


2.73 


4.86 


11.32 


160 


3.06 


6.55 


17.85 


2.34 


4.19 


9.22 


165 


2.63 


5.68 


14.66 


2.01 


3.63 


7.57 


170 


2.28 


4.95 


12.03 


1.73 


3.12 


6.26 


175 


1.96 


4.29 


10.03 


1.49 


2.71 


5.13 


180 


1.70 


3.74 


8.24 


1.28 


2.34 


4.23 


185 


1.47 


3.26 


6.86 


1.11 


2.04 


3.51 


190 


1.28 


2.85 


5.65 


0.96 


1.77 


2.88 


195 


1.11 


2.48 


4.74 


0.83 


1.54 


2.38 


200 


0.96 


2.17 


3.89 


0.72 


1.34 


1.98 




Q 2 = m 2 /4 fi™ = 0.067m, /i^ s = 0.13m 


150 


7.10 


13.16 


12.44 


5.41 


8.84 


7.42 


155 


6.13 


11.44 


10.21 


4.65 


7.63 


6.04 


160 


5.30 


9.91 


8.35 


3.99 


6.58 


4.91 


165 


4.58 


8.65 


6.86 


3.44 


5.69 


4.03 


170 


3.95 


7.49 


5.64 


2.97 


4.94 


3.28 


175 


3.44 


6.57 


4.65 


2.57 


4.29 


2.71 


180 


2.98 


5.72 


3.84 


2.23 


3.75 


2.25 


185 


2.60 


5.02 


3.20 


1.93 


3.25 


1.84 


190 


2.27 


4.40 


2.65 


1.68 


2.84 


1.52 


195 


1.98 


3.85 


2.19 


1.46 


2.48 


1.25 


200 


1.73 


3.38 


1.82 


1.28 


2.17 


1.04 



Table 2: The hadronic it production cross sections in pp collisions at v^S* = 
1.8 TeV for s cut = 5s 4imin and 9* = 90°. The effective /i from |2| is given for 
each Q 2 . 
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NLL 


LL 


m (GeV/c 2 ) 


^ IS (Pb) 


^ S (Pb) 


°99 (Pb) 


^ IS (Pb) 


< S (Pb) 


°99 (Pb) 




Q 2 = m 2 /^ IS = 0.19m, /i^ s = 0.33m 


150 


4.40 


6.30 


1.61 


3.61 


4.90 


1.24 


155 


3.78 


5.43 


1.29 


3.10 


4.23 


0.98 


160 


3.25 


4.69 


1.03 


2.66 


3.64 


0.78 


165 


2.80 


4.06 


0.83 


2.28 


3.13 


0.63 


170 


2.42 


3.53 


0.67 


1.97 


2.71 


0.50 


175 


2.10 


3.06 


0.54 


1.70 


2.34 


0.40 


180 


1.82 


2.67 


0.44 


1.47 


2.03 


0.33 


185 


1.58 


2.32 


0.35 


1.27 


1.76 


0.26 


190 


1.37 


2.02 


0.29 


1.10 


1.53 


0.21 


195 


1.19 


1.76 


0.23 


0.96 


1.33 


0.17 


200 


1.04 


1.55 


0.19 


0.83 


1.16 


0.14 




Q 2 = 4m 2 /iP IS = 0.27m, /i^ s = 0.52m 


150 


3.55 


5.36 


1.70 


2.91 


4.11 


1.24 


155 


3.05 


4.61 


1.36 


2.48 


3.52 


0.99 


160 


2.62 


3.98 


1.09 


2.13 


3.03 


0.79 


165 


2.25 


3.44 


0.87 


1.83 


2.61 


0.63 


170 


1.94 


2.97 


0.71 


1.57 


2.25 


0.51 


175 


1.68 


2.58 


0.57 


1.35 


1.94 


0.41 


180 


1.45 


2.24 


0.46 


1.17 


1.68 


0.33 


185 


1.26 


1.95 


0.38 


1.01 


1.46 


0.27 


190 


1.09 


1.69 


0.31 


0.87 


1.26 


0.22 


195 


0.95 


1.48 


0.25 


0.75 


1.09 


0.18 


200 


0.82 


1.28 


0.20 


0.69 


0.96 


0.14 




Q2 = m 2/ 4 ^ ^dis = o.l3m, /i^ s = 0.21m 


150 


5.27 


7.10 


1.57 


4.40 


5.68 


1.25 


155 


4.54 


6.14 


1.25 


3.77 


4.89 


1.00 


160 


3.92 


5.31 


0.99 


3.25 


4.23 


0.80 


165 


3.39 


4.62 


0.80 


2.80 


3.65 


0.64 


170 


2.94 


4.02 


0.64 


2.41 


3.16 


0.51 


175 


2.55 


3.49 


0.52 


2.08 


2.74 


0.41 


180 


2.22 


3.05 


0.42 


1.81 


2.38 


0.33 


185 


1.93 


2.67 


0.34 


1.57 


2.07 


0.27 


190 


1.68 


2.33 


0.28 


1.37 


1.80 


0.22 


195 


1.47 


2.04 


0.23 


1.18 


1.57 


0.18 


200 


1.28 


1.78 


0.18 


1.03 


1.37 


0.14 



Table 3: The hadronic tt production cross sections in pp collisions at v^S* = 
1.8 TeV for s cut = 20s 4jmin and 9* = 90°. The effective /x from M is given 
for each Q 2 . 
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NLL 


LL 


Pi a h (GeV/c) 


^ IS (nb) 


</ (nb) 


agg (nb) 


o%* (nb) 


<f (nb) 


(Tgg (nb) 




Q2 = m 2 ^ ^ms = o.37m, fi™ s = 0.51m 


200 


0.0049 


0.0081 


0.034 


0.0021 


0.0031 


0.014 


300 


0.049 


0.091 


0.410 


0.023 


0.037 


0.172 


400 


0.169 


0.326 


1.54 


0.083 


0.140 


0.658 


500 


0.373 


0.734 


3.58 


0.190 


0.325 


1.56 


600 


0.653 


1.30 


6.52 


0.344 


0.591 


2.90 


700 


1.00 


2.01 


10.44 


0.537 


0.927 


4.64 


800 


1.40 


2.84 


15.00 


0.764 


1.32 


6.76 


900 


1.84 


3.72 


20.23 


1.02 


1.77 


9.20 


1000 


2.31 


4.70 


26.18 


1.30 


2.25 


11.86 




Q 2 = Am 2 =>- /i° IS = 0.52m, /ijf s = 0.82m 


200 


0.0036 


0.0079 


0.041 


0.0016 


0.0029 


0.015 


300 


0.036 


0.082 


0.482 


0.017 


0.032 


0.183 


400 


0.124 


0.286 


1.81 


0.060 


0.118 


0.702 


500 


0.273 


0.636 


4.32 


0.138 


0.269 


1.69 


600 


0.479 


1.12 


8.01 


0.248 


0.484 


3.20 


700 


0.732 


1.72 


12.87 


0.387 


0.756 


5.15 


800 


1.02 


2.41 


18.86 


0.551 


1.07 


7.60 


900 


1.35 


3.17 


25.91 


0.735 


1.43 


10.42 


1000 


1.70 


4.01 


33.68 


0.941 


1.83 


13.69 




Q2 = m 2/ 4 ^ ^dis = o.26m, /i^ s = 0.32m 


200 


0.0041 


0.0054 


0.026 


0.0019 


0.0022 


0.012 


300 


0.050 


0.076 


0.354 


0.024 


0.033 


0.167 


400 


0.188 


0.300 


1.37 


0.096 


0.139 


0.665 


500 


0.433 


0.713 


3.20 


0.232 


0.341 


1.59 


600 


0.783 


1.30 


5.78 


0.429 


0.640 


2.94 


700 


1.21 


2.04 


9.11 


0.682 


1.03 


4.69 


800 


1.72 


2.92 


12.96 


0.985 


1.49 


6.79 


900 


2.27 


3.90 


17.37 


1.33 


2.02 


9.16 


1000 


2.87 


4.95 


22.29 


1.70 


2.60 


11.81 



Table 4: The hadronic bb production cross sections in pp interactions for 
s cut = 1.4s 4 ,min, m = 4.75 GeV/c 2 and 9* = 90°. 
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NLL 


LL 


Pi a h (GeV/c) 


^ IS (nb) 


</ (nb) 


agg (nb) 


^ IS (nb) 


<f (nb) 


<7gg (nb) 




Q2 = m 2 ^ ^ms = o.34m, fi™ s = 0.49m 


200 


0.0094 


0.020 


0.248 


0.0029 


0.0049 


0.051 


300 


0.087 


0.206 


2.75 


0.029 


0.055 


0.594 


400 


0.287 


0.700 


9.82 


0.103 


0.199 


2.19 


500 


0.617 


1.54 


22.43 


0.232 


0.452 


5.06 


600 


1.06 


2.67 


40.53 


0.413 


0.811 


9.21 


700 


1.59 


4.05 


63.18 


0.639 


1.26 


14.67 


800 


2.20 


5.62 


89.85 


0.902 


1.78 


20.84 


900 


2.85 


7.32 


120.6 


1.20 


2.37 


28.00 


1000 


3.56 


9.18 


154.2 


1.51 


2.99 


36.60 




g 2 = 4m 2 =>- /i° IS = 0.48m, /i^ s = 0.77m 


200 


0.0063 


0.017 


0.301 


0.0019 


0.0042 


0.057 


300 


0.059 


0.171 


3.41 


0.020 


0.045 


0.665 


400 


0.193 


0.575 


12.56 


0.071 


0.159 


2.51 


500 


0.416 


1.25 


29.35 


0.159 


0.356 


5.96 


600 


0.716 


2.16 


54.00 


0.284 


0.635 


11.10 


700 


1.08 


3.30 


85.96 


0.442 


0.985 


17.73 


800 


1.49 


4.58 


125.7 


0.626 


1.39 


26.18 


900 


1.95 


5.97 


170.1 


0.829 


1.84 


35.82 


1000 


2.45 


7.47 


220.3 


1.05 


2.34 


46.82 




g 2 = m 2 /4 fj$ ls = 0.24m, /ijf 3 = 0.31m 


200 


0.0099 


0.016 


0.201 


0.0030 


0.0041 


0.046 


300 


0.106 


0.200 


2.35 


0.035 


0.056 


0.571 


400 


0.369 


0.735 


8.33 


0.132 


0.218 


2.10 


500 


0.813 


1.66 


18.63 


0.306 


0.516 


4.80 


600 


1.41 


2.94 


32.98 


0.557 


0.948 


8.63 


700 


2.14 


4.51 


50.25 


0.868 


1.49 


13.36 


800 


2.97 


6.33 


70.60 


1.24 


2.13 


19.03 


900 


3.85 


8.26 


92.73 


1.65 


2.85 


25.16 


1000 


4.78 


10.32 


116.5 


2.10 


3.63 


32.12 



Table 5: The hadronic bb production cross sections in pp interactions for 
s cut = 1.2s 4 ,min, m = 4.75 GeV/c 2 and 9* = 90°. 
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NLL 


LL 


Pi a h (GeV/c) 


^ IS (nb) 


</ (nb) 


agg (nb) 


^ IS (nb) 


<f (nb) 


<jgg (nb) 




Q2 = m 2 ^ ^ms = o.39m, fi™ s = 0.54m 


200 


0.0035 


0.0050 


0.015 


0.0017 


0.0023 


0.0071 


300 


0.037 


0.061 


0.187 


0.019 


0.029 


0.094 


400 


0.131 


0.225 


0.723 


0.072 


0.111 


0.373 


500 


0.295 


0.517 


1.72 


0.166 


0.262 


0.910 


600 


0.525 


0.932 


3.20 


0.303 


0.482 


1.70 


700 


0.811 


1.46 


5.14 


0.478 


0.765 


2.76 


800 


1.14 


2.06 


7.46 


0.685 


1.10 


4.05 


900 


1.51 


2.74 


10.18 


0.918 


1.48 


5.54 


1000 


1.92 


3.48 


13.20 


1.17 


1.90 


7.24 




Q 2 = Am 2 =>- /i° IS = 0.56m, /ijf 3 = 0.85m 


200 


0.0028 


0.0054 


0.018 


0.0013 


0.0023 


0.0079 


300 


0.029 


0.059 


0.216 


0.015 


0.027 


0.099 


400 


0.102 


0.209 


0.825 


0.054 


0.097 


0.388 


500 


0.228 


0.471 


1.99 


0.126 


0.227 


0.945 


600 


0.402 


0.837 


3.72 


0.227 


0.409 


1.78 


700 


0.620 


1.29 


6.02 


0.356 


0.644 


2.93 


800 


0.874 


1.83 


8.85 


0.509 


0.920 


4.31 


900 


1.16 


2.43 


12.17 


0.684 


1.23 


5.93 


1000 


1.46 


3.06 


15.94 


0.873 


1.58 


7.88 




g 2 = m 2 /4 fj$ ls = 0.28m, /ijf 3 = 0.34m 


200 


0.0025 


0.0029 


0.010 


0.0013 


0.0014 


0.0057 


300 


0.033 


0.045 


0.159 


0.018 


0.023 


0.090 


400 


0.133 


0.188 


0.646 


0.076 


0.101 


0.376 


500 


0.315 


0.461 


1.57 


0.188 


0.255 


0.940 


600 


0.580 


0.866 


2.93 


0.356 


0.489 


1.78 


700 


0.915 


1.38 


4.69 


0.576 


0.801 


2.89 


800 


1.31 


1.99 


6.81 


0.837 


1.17 


4.25 


900 


1.76 


2.70 


9.26 


1.14 


1.60 


5.85 


1000 


2.24 


3.46 


11.92 


1.47 


2.08 


7.58 



Table 6: The hadronic bb production cross sections in pp interactions for 
s cut = 1.6s 4 ,min, m = 4.75 GeV/c 2 and 9* = 90°. 
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